Polymer drift in a solvent by force acting on one polymer end 
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We investigate the effect of hydrodynamic interactions on the non-equilibrium drift dynamics of 
an ideal flexible polymer pulled by a constant force applied at one end of the polymer using the 
perturbation theory and the renormalization group method. For moderate force, if the polymer 
elongation is small, the hydrodynamic interactions are not screened and the velocity and the lon- 
gitudinal elongation of the polymer are computed using the renormalization group method. Both 
the velocity and elongation are nonlinear functions of the driving force in this regime. For large 
elongation we found two regimes. For large force but finite chain length L the hydrodynamic inter- 
actions are screened. For large chain lengths and a finite force the hydrodynamic interactions are 
only partially screened, which in three dimensions results in unusual logarithmic corrections to the 
velocity and the longitudinal elongation. 

PACS numbers: 05.40.-a, 36.20.-r, 61.41.+e 



I. INTRODUCTION 

Understanding the dynamics of polymers in solutions is 
of great importance in many bio-chemical and industrial 
process. In dilute solution, hydrodynamic interaction be- 
tween different parts of the polymer chain plays a crucial 
role in the determination of dynamic properties In 
the classical theories which go back to Kirkwood and 
Zimm 0, Q , the main focus of research was to study the 
quantities such as the diffusion constant, the dynamic in- 
trinsic viscosity, and the dynamic structure factor of the 
polymer. Recent advances in experimental techniques 
make it possible to explore the motion of individual poly- 
mer under hydrodynamic flow, thermal noise, or external 
fields. For instance, Chu et al. 0, H, [f| experimentally 
studied the behavior of polymers in different flow condi- 
tions with the emphasis on biological applications. See 
also related theoretical studies 

In this paper we study drift of an ideal flexible polymer 
pulled by a constant driving force applied at one poly- 
mer end. Classical theory of dragging objects dates back 
to Einstein's Brownian motion and Stokes' theory in the 
hydrodynamic medium, which yields the velocity of the 
dragged object, v B = F/fo with / = Ci] s R. Unlike the 
drag of these simple objects, the dynamics of the poly- 
mer pulled at polymer end is more complicated due to 
the interplay between the internal degrees of freedom of 
polymer, hydrodynamic interactions, and inhomogeneity 
of chain configurations, which is due to the pulling force 
acting at polymer end. The drift by a force applied at 
polymer end is different from the case, if the external 
force is homogeneous and acts on all monomers of the 
chain. 

In neglecting the hydrodynamic interactions i.e. in the 
Rouse theory the average velocity and the chain elonga- 
tion are given by expressions 



(r z (0,t)-r*(L,t)) 



FLl 
2dk B T' 



(2) 



where d is the dimensionality, F is the pulling force, and 
/o denotes the monomer friction coefficient. For small 
forces the polymer is expected to have the shape of coil 
with the consequence that the hydrodynamic interaction 
is not screened (Zimm regime), and the velocity is de- 
scribed by Stokes formula 
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with Rh ~ WN being the hydrodynamic radius of the 
polymer, and rj s is the solvent viscosity. The Zimm be- 
havior ([3]) is described in the framework of the renormal- 
ization group method as due to the renormalization of the 
friction coefficient, which becomes scale dependent in the 
regime of strong hydrodynamic interaction (lfj|-p^|. 

For moderate forces a crossover between ([TJ and ([3]) 
as well as a nonlinear dependence on the pulling force 
is expected to apply. For large forces the effect of the 
hydrodynamic interaction is expected to be weak due to 
polymer stretching, so that a Rouse-type behavior should 
be observed. 

In the present work we calculate the first-order correc- 
tions to the velocity and the longitudinal elongation. The 
renormalization group method allows us to establish the 
behavior of the quantities under consideration for mod- 
erate forces. We also show that in d = 3 the first-order 
corrections to the velocity and the longitudinal size of the 
polymer in powers of the hydrodynamic interactions loga- 
rithmically depend on the parameter f3Fy^lL/2d, which 
demonstrates that hydrodynamic interactions are only 
partially screened for long polymer pulled at finite force. 

The paper is organized as follows. After an introduc- 
tion to the formalism in section[TIJ we perform the pertur- 
bational computation of the average velocity in section 
Mil and the longitudinal polymer size under the drift in 
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section HVl In section [V] the velocity and the longitudi- 
nal size are analyzed for small and large force using the 
renormalization group analysis. The work is summarized 
in section [VTl 



II. THE FORMALISM 

The drift dynamics of a polymer chain is described by 
the Kirkwood diffusion equation, which in the continuum 
limit is given by [l[ 



L A f\ 6 

CIS2 / dS\-^-, r x 



dP({r(s)},t) 

dt J J ~'5r(s 1 

[D 5^6( Sl - s 2 ) + T""(r(s 1 ) - r(s 2 ))] x 



5r(s 2 ) 



Sr(s 2 ) 



P. 



(4) 



The Oseen tensor (multiplied by ksT) in d-dimensions is 
defined through its Fourier transform as 

T^(r( Sl ) - r(s 2 )) = 

fcsT f d d q 1 g^g" ici M sl )-r(s 2 )) 

Vs J (2^v l q i >e 

The free energy (multiplied by fcgT) and the pulling force 
applied in z-direction are given by 



, F^ t (v(s))=S fJ ' z Fd(s). (5) 



The Rouse modes £ fc are the Fourier coefficients in the 
expansion of the position vector r(s) (0 < s < L) accord- 

oo 

ing to r(s) = ^2 Qsk^k where the basis functions 

fc=0 
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k = 

l i ' k = 1,2,3,- ■■ 



2 cos 7r ' M 
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are orthogonal and complete in the interval (0, L). 

The Kirkwood diffusion equation can be rewritten in 
terms of the normal coordinates 



dP(Z,t;£°,t ) 
dt 



= UP + UP, 



(7) 



where £ is abbreviation for , and the operators U and 
Li are given by 

U = D Vl + AAfcV^ - D pF^Q ok Vl (8) 



Li = (V^ + A m C) - K f ko^ (9) 

with abbreviations V&=£/££;fc, 

T fem = / <&2 / dsiQ^feT^^Si) - r(s 2 ))Q S2m , 
Jo Jo 

f% = [ d Sl Q Slk T^(r( Sl )-r(0)), 



and r"(*i) - r^(s 2 ) = £ (Q Sl „ - Q S2 „)££. 

n=0 

In the following we consider the hydrodynamic inter- 
action as perturbation, and use the Rouse chain pulled 
by a force at one end as a reference state. The solution 
of Eq. ([7]) with Li = gives the transition probability 
density of the reference Rouse state 
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(10) 



with 



a fe (t) = exp(-L> A fe i), /fc(i) 



1 - cxp(-2£>oA fc t) 
2A fc 



where Z?o = kgTl/ /o is the monomer self-diffusion coef- 
ficient, and A fe = j (^) are the Rouse eigenvalues. 

To perform perturbational computations of the veloc- 
ity, the longitudinal size, etc. of the polymer chain un- 
der the pulling force, we rewrite the Kirkwood diffusion 
equation as an integral equation 



P(U;f,t )=P ^,f,f,to)+ / dt'x 

Jt 

DZ'Poit t; C t')U{C)P{e, t'; e, to). (11) 



Iteration of (|f I p results in the perturbation expansions in 
powers of the hydrodynamic interaction, which is sym- 
bolically written: P = P + PqL.P = P + PqL^o + 
PoUiPoLiP = ■ ■ ■ ■ 



III. PERTURBATIONAL COMPUTATION OF 
THE VELOCITY 

Our aim is now to calculate the velocity of the polymer 
chain up to the first order in hydrodynamic interaction 
using Eq. (fTTj) . The center of mass r c of the polymer 
chain is expressed through the center of mass Rouse mode 
£fe=o according to r c = £ k=0 /VZ = (l/L)J Q L dsr(s). 
Hence, the computation of the average velocity reduces 
to the computation of the expectation value of the zero 
mode (r c (i)) = (v c (t)) t: 



(*k=o(«)) 



J D^ k=0 p{ct-e,h). 



(12) 



The computation of (TT2")) to the zero order in powers of 
the hydrodynamic interaction (i.e. in the Rouse model) 
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using the unperturbed distribution function (|10[) yields 
in the steady state 



(S=o(*)> ( 



IF 

7vz 



t, 



where we take the force, which is applied at s = 0, to 
be directed along the z-axes, F — F z . For the velocity 
we obtain the expression (v*(t)) = F/(fN), which is in 
agreement with Rouse theory. The first-order correction 
is derived from Eq. (jTTJ) by replacing P on the right-hand 
side by Pq as 



(a=o(*)>!= fdf r d S2 r d S1 1 ,, " q 

Jt JO JO 



T zz (q)Q sl o I Dd'e^So^' 



PFQoo x 



Pott'MSo - $ ~ D f3FQ 00 (t' - i )) 



(13) 



where T^{q) 
over £' yields 



k&L^^ _ 2_g_). The integrations 
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Q sin Q S2n 
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The evaluation of the sums in Eq. (fT4")l gives 

7 1zz (g)Q Sl0 /3^Qooexp(^ J F 1 z /?6( Sl , S2 )-a g 2 ) (15) 
where the quantities a and b are defined as follows 

(Qsin Q S2ft) ^ 

,r J 

Qsxn Qs 2 n ^ ' / « /-, s l + s 2 \ 
7 Qon = -j(s 2 - Si)(l - 



° = ^ 2A„ 

Qsin Q S2T1 
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Performing the straightforward integration over q we fi- 
nally obtain 



1A 
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where the function A(y) and its argument y are respec- 
tively given by 



A{y) 



(d-l)(d-2) (d-l)(rf-2-2y 2 



4-d /d-4 



d- 2 



(17) 



(18) 



T(a,z) = J^°dtt a ~ 1 e~ t is the incomplete gamma func- 
tion. The function A(y) behaves for small and large ar- 
guments as 



d(d-2) ~ d(d+2)V + {d+2)(d+i)y H ' y^- 1 -^ 



4(^-1) 



(d - 1)( y +2 - d) r(^) + ^e^ 2 +..., y»l. 

(19) 

Integrations over x\ and x 2 in Eq. (|16p . which can be 
carried out in the limit F — > 0, yield 



A(0) 



(6-d)(4-d) 



A(0) 



d~4 



4-d 



A(0). (20) 



Note that the behavior of the first-order correction to 
the velocity in the vicinity of four dimensions, which is 
given by the last expression, plays an important role in 
the renormalization group analysis of the velocity. 



IV. PERTURB ATIONAL COMPUTATION OF 
THE LONGITUDINAL SIZE OF THE POLYMER 



In this section we calculate the longitudinal size of the 
polymer up to the first order in hydrodynamic interac- 
tion. The longitudinal size of the polymer is expressed 
through the Rouse modes according to 



(0,t)-r*(L,t)} = J-'£(l-(-l) k ){e k (t)}. (21) 



k=l 



The computation to the zero order in hydrodynamic in- 
teraction (i.e. in the Rouse model) in the steady state 
yields 



(r z (0,i)-r*(i,i)) o 



FLl 
2dk B T' 
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The first-order correction to (£|(t)) is obtained using 
Eq. QIJ as 

(£ s z k (t)) 1 = (3F I dt'a k (t-t') ds 2 d Sl x 

Jt Jo JO 

d d q 

J^yiQs 2 kT zz (q)Q Sim Qom |m=0 X 



exp \iq F(3 



A, 



I On- 



E 



2 A n 



(22) 



The integration over q and summations over k and n are 
performed similar to those in the preceding section. The 
longitudinal elongation reads 



(r z (0,t)-r z (L,t)) 



ii (1 _IA x 
2dk B T v 2 2 drj s 



d/2 



2,1) L2 - d/2 ^Fff d ) + --) 



(23) 



with 



B r ((3Fy/lL/2d) = dx 2 dx x 

JO Jo \ x 2 — x l) ' 



(2x 2 - l)A(y) 



The computation of integrals in Eq. 
for a small force yields 

£ r (F->0) 



2-1 

(24) 
over xi and x-i 



(8-d)(6-cZ) 



A(0). 



The finiteness of the latter at four dimensions means in 
the context of the renormalization group method that 
the parameter F, L, I, and T appearing in the prefactor 
of Eq. ([23| do not renormalize. The only quantity which 
renormalize is the monomer friction coefficient. The hy- 
drodynamic interaction results in decreasing the size of 
the polymer. 



V. RESULTS 

Our aim is now to study the behavior of the veloc- 
ity and the longitudinal size of the polymer chain as a 
function of the driving force F and the chain length L. 



A. Small elongation 

The first-order perturbational correction to the veloc- 
ity is the starting point to perform the renormalization 
group (RG) analysis, which enables one to take into ac- 
count the effect of hydrodynamic interaction beyond the 
first-order. The basic observation to apply RG is that 
the integral in (fT6|) diverges logarithmically in four di- 
mensions (i.e. the critical dimension is four; for d > 4 the 



hydrodynamic interactions become irrelevant). This di- 
vergence manifests itself in d < 4 dimensions as 1/(4— d) 
pole. To regularize the theory these poles in perturba- 
tion expansions have to be removed by an appropriate 
renormalization of the friction coefficient. In the limit of 
a small pulling force the renormalizcd friction coefficient 
is derived from (fTHl) as 



(25) 



/ = /o(l-^o(i e/2 -A^) + ...) 
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where £o = (f /(r] s d))(d/(2Trl)) d / 2 is the expansion pa- 
rameter of perturbation expansions in powers of the hy- 
drodynamic interactions, and e = 4 — d. The ultraviolet 
cutoff A in ([25| is introduced to enable the limit to four 
dimensions. The cutoff excludes the hydrodynamic in- 
teractions between monomers separated along the chain 
by the contour length less than A, |si — S2 1 < A. The 
renormalization of the friction coefficient ()25|1 obtained 
from Eq. (|16p coincides with that obtained by studying 
different problems in polymer dynamics [l(| -P3- Note 
that in the absence of the excluded volume interaction, 
only friction coefficient renormalizes. 

In addition to the renormalization of the friction coeffi- 
cient one should consider the renormalization of the cou- 
pling constant (i.e. expansion parameter) controlling the 
strength of the hydrodynamic interaction. The inspec- 
tion of Eq. (|16jl or Eq. |23|) yields the bare dimensionless 
expansion parameter as 



-o = 4(/y) d/2i£/2 - 

rj s d 27rt 



(26) 



It follows from Eq. (|26|) that the coupling constant renor- 
malizes in the same way as the friction coefficient. 

The renormalization group is based on the observa- 
tion that the regularization, i.e. the elimination of the 
1 /e-poles from perturbation expansions can be performed 
step by step by changing the cutoff A — * A' — * A" — > 
■ • • — * A m . The renormalization of the friction coefficient 
and the strength of the hydrodynamic interaction due to 
an infinitesimal change of the cutoff are given by differ- 
ential equations. To the one-loop order one obtains 



v df 3 

dW = 4 W ' 



, dw e 3 2 

A av^2 w "I w 



(27) 



(28) 



with the dimensionless effective coupling constant w = 
(f/(rj s d))(d/(27rl)) d / 2 X' E / 2 . The solutions of Eqs. $7$ 
and (2SI read 



/o 



l + iCo(AL /2 -A^)' 



w = -?-{— ) d/2 X e/2 . 
Tj s d 2irl m 



(29) 



(30) 
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It follows from Eq. ([29)1 and (|30|) that at large A m , w 
approaches the fixed-point value w* — 2e/3, which corre- 
sponds to the zero point of the Gell-Mann - Low function 
/3(w). At the fixed-point the effective friction coefficient 
depends on X m as power law A m e ^ 2 . At low forces A m 
is equal to L, so that the renormalized (effective) fric- 
tion coefficient scales as / = w*rj s 

jd/2£- e /2_ The drift 

velocity behaves consequently as 



F 



F 



(31) 



which agrees in d = 3 with the Stokes formula ([3]). 

The inspection of the first-order correction to the lon- 
gitudinal size of the polymer yields that it is finite in four 
dimensions. This is what is expected, because the fric- 
tion coefficient, which is the only quantity to renormalize, 
does not appear in the zero-order correction to the lon- 
gitudinal size of the polymer. Thus, the RG prediction 
for the longitudinal size consists in replacing the bare ex- 
pansion parameter in Eq. (|23|) by the renormalized one. 
According to Eq. (|3"0|) the first-order correction is of order 
e and is thus small in the vicinity of four dimensions. 

For small forces the velocity and the elongation are 
given in the renormalized theory by expressions 



W {1 



-wB v (l3F^/lL/2d) 



(32) 



(r z (0,t)~r z (L,t)) = 



FLl 



1 



The function B v (y) in 



JWfcflT ( 1 -4^iVtW)+-). 

(33) 

is defined by B v (y) = 
The effective expansion pa- 
O(e)), so that the ex- 



(x 2 -X 1 ) d / 2 - 1 ■ 

rameter w is a small number 
pansions (|32[) and (f33|) are reliable. The expansion of the 
functions B v (z) and B r (z) for small z are given by 



B v (z) = 



1136 



z- + 



14175 

64 464 
B r (z) = — - 



19366 

7882875 

10786 



(34) 

-z 4 + • • • . (35) 
45 31185 19144125 v ; 

Expressions (|32j|33p demonstrate that the velocity and 
the longitudinal polymer size show the non-linear re- 
sponse i.e. non Stokes behavior at moderate pulling 
forces. 

Note that the Zimm result for the velocity Q can be 
also obtained in a simple way following the Kirkwood 
type approach to the diffusion constant [2} ■ That consists 
in neglecting one in the bracket of (fl6|) . The friction 
coefficient then drops, and the dependence on N is the 
same as ©. 

B. Large elongation 

We now will study the first-order corrections to the 
velocity and the longitudinal elongation for large argu- 
ment (3Fy / lL/2d. To estimate the integrals in (JTJ 



for large forces we use the asymptotic expression of A(y) 
given by Eq. (|19[) . To prevent the divergence in the inte- 
gral over t = X2 — x\ from the integration region t = 0, 
we introduce an cutoff to — 6/Llf3 2 F 2 . The evaluation of 
the integrals over t and can be performed analytically 
and gives for d — 3 

with a numerical constant c v . We would like to stress 
that the integral in Eq. (|16p computed numerically can 
be fitted with (In?/) 1 - 75 /y instead of the estimate given 
by Eq. (f3"6"]) . The difference is due to the complicated 
form of the expression for y in Eq. (fT5| . 

For large F and finite L we arrive at the Rouse result. 
However, for large L and finite F the correction increases 
logarithmically with L, and will become large for large 
L. Unfortunately, there are no analytical means to study 
the effect of the whole perturbation expansion on v z in 
this case. The extrapolation of (f36|) a la Kirkwood, i.e. 
the disregard of one in the bracket of Eq. (|31)|) , yields 



(37) 



Eq. (|37|) shows that in this regime (finite F and large L) 
the hydrodynamic interactions determine the behavior of 
the polymer. The friction coefficient drops in the expres- 
sion of the velocity in this regime. We expect that the 
hydrodynamic interactions remain important due to the 
thermal fluctuations of the pulled chain. The presence of 
ksT in (|3T|) supports this statement. Note that the loga- 
rithmic dependencies of v z c on the force does not allow to 
write Eq. (|37|) in the form of Eq. HJ with some effective 
friction coefficient. Eq. (|37[) can be formally obtained 
from Eq. |T]) by simultaneous replacements 



fo^ f ~l Vs /\n((3F(Ll) 1/2 ), F 



ksT (3F{Ll) 1/2 

m . 

I 3 

(38) 

The "renormalized" friction coefficient in Eq. (|3"5)) resem- 
bles to some extent the friction coefficient, Lr\ s j ln(L/Z), 
of a slender body (a cylinder of length L and cross section 
radius I, L ^> I) in a flow UGl (see also U). However, 
the elongated polymer chain is different from a rod, and 
the comparison is only qualitativ. 

A similar computation of the first-order correction to 
the chain size yields 



FLl 



(r z (Q,t)-r z (L,t)) = 
X V ' ' V ' " 2dk B T 

(l n £M^_ 2 )ln (f3F(Ll) 1/2 ) 



(l-C^oL 1 ' 2 : 



(3F (LI) 



1/2 



(39) 



Note that the integral in Eq. (f2"4")) computed numerically 
can be fitted with (\ny) 2Ar /y instead of the analytical 
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estimate in Eq. (|39|) . For large F and finite L, we ar- 
rive at the Rouse result. On the contrary, for large L 
and finite F the first-order correction to the longitudi- 
nal size increases logarithmically with L. This results in 
a decrease of the longitudinal size of the polymer. It is 
difficult, due to the absence of an analytical method to 
make predictions on the total effect of the hydrodynamic 
interaction for polymer elongation. Note that the cor- 
rections in Eqs. (fl"6l [23)) at F 7^ imply the nonlinear 
response of the polymer on the pulling force. 

VI. SUMMARY 

To summarize, we have studied the drift of an ideal 
polymer driven by a constant force applied at one poly- 
mer end using perturbation expansions in powers of hy- 
drodynamic interactions. For moderate force, where 
the hydrodynamic interactions are not screened and the 
polymer elongation is small, the renormalization group 
method permits to compute the velocity and the longi- 
tudinal elongation of the polymer. These quantities are 
nonlinear functions of the driving force in this regime. 
For large elongation we found two regimes. For large 



force but finite chain length L the hydrodynamic inter- 
actions are screened, so that both the velocity and the 
longitudinal polymer elongation are given by the corre- 
sponding results of the Rouse theory. For large L but 
finite force, the regime which we have studied for d = 3, 
the hydrodynamic interactions are partially screened. 
The first-order corrections to v z and (r z (0,t) — r z (L,t)}, 
increase logarithmically with L. Following the Kirk- 
wood's treatment of the diffusion constant we make a pre- 
diction for the velocity, Eq. (j3"o]l . beyond the first-order 
of the perturbation theory. It would be of great inter- 
est to check experimentally these predictions by pulling 
polymer in solvent using optical tweezers. 
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